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Abstract 

The four most stable structures of Ni^r clusters with N from 2 to 150 have been determined 
using a combination of the embedded-atom method in the version of Daw, Baskes and Foiles, the 
variable metric/ quasi- Newton method, and our own Aufbau/Abbau method. A systematic study of 
energetics, structure, growth, and stability of also larger clusters has been carried through without 
more or less severe assumptions on the initial geometries in the structure optimization, on the 
symmetry, or on bond lengths. It is shown that cluster growth is predominantly icosahedral with 
islands of fee, tetrahedral and decahedral growth. For the first time in unbiased computations 
it is found that M147 is the multilayer (third Mackay) icosahedron. Further, we point to an 
enhanced ability of fee clusters to compete with the icosahedral and decahedral structures in the 
vicinity of iV = 79. In addition, it is shown that conversion from the /icp/anti-Mackay kind of 
icosahedral growth to the /cc/Mackay one occurs within a transition layer including several cluster 
sizes. Moreover, we present and apply different analytical tools in studying structural and energetic 
properties of such a large class of clusters. These include means for identifying the overall shape, the 
occurrence of atomic shells, the similarity of the clusters with, e.g., fragments of the fee crystal or of 
a large icosahedral cluster, and a way of analysing whether the iV-atom cluster can be considered 
constructed from the (N — l)-atom one by adding an extra atom. In addition, we compare in 
detail with results from chemical-probe experiment. Maybe the most central result is that first for 
clusters with N above 80 general trends can be identified. 
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I. INTRODUCTION 



Clusters are important materials both from the point of view of basic research and from an 
application point of view. Their, partly controllable, unique physical and chemical properties 
can be related to the large surface-to-volume ratio or to finite-size or quantum-confinement 
effects. These properties make them interesting for use in, e.g., nanoelectronics and in 
catalysis. 

The central issue is to understand and predict how the properties of interest depend on the 
size of the cluster. However, addressing this question is complicated by a number of serious 
problems. First, clusters contain typically from some 10s to several 100s of atoms, so that 
they are neither small, well-characterized molecules, nor macroscopic, approximately infinite 
crystals. Moreover, in experimental studies, one often considers an ensemble of more or less 
mono-disperse clusters whose size distribution is largely unknown. Second, experimental 
and theoretical studies focus often on different systems: in experiment the clusters may be 
dispersed in some solvent and may possess surfactants and only sometimes are produced in 
the gas phase, whereas isolated clusters in this gas phase almost exclusively are the subject of 
theoretical studies. Third, in theoretical works the fairly large size of the individual clusters 
makes it necessary to invoke one or more approximations: either highly accurate methods 
on selected clusters with pre-chosen structures are applied, or less accurate, parameterized 
methods are used on more different clusters, but also then the structures are often chosen 
'reasonably'. 

In this work we shall present results of our theoretical studies on Ni^y clusters with iV 
up to 150. Our approach is based on the approximate embedded-atom method (EAM) 
in calculating the total energy of a given structure and we shall show that this method 
provides a good compromise between accuracy and computational speed. Moreover, we 
have used our own Aufbau/Abbau method in determining the structures of the lowest total 
energies. Ni clusters provide one of the most studies sets of clusters and, therefore, is an 
excellent system for exploring new theoretical or experimental approaches. This is our main 
reason for focusing on that system. Moreover, through our unbiased structure-optimization 
approach we are able to identify the structures of more low-energy structures for each size, 
and, simultaneously, obtain structure and total energy for all values of N up to iV = 150, 
whereby size-specific properties can be extracted. 
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Once the structure of a given cluster has been chosen, the variational principle allows 
for a systematic improvement of the calculated electronic properties. A similar approach 
does not exist for the structural properties where one, instead, has to use more or less 
biased methods for determining the structure of the lowest total energy. Here, three issues 
make such theoretical studies very demanding: first, the computational efforts for a single 
structure scales with N to some power from 2 and upwards; second, the number of structural 
parameters that have to be determined scale linearly with N; and, third, the number of 
meta-stable structures scales faster than any power of iV.— Therefore, any theoretical study 
employs one or more approximation that may or may not be crucial. 

Accordingly, parameter-free calculations on Nijv cluster o 2 ' 3 ' 4 ' 5 ' 6 ' 7 have been applied mainly 
for smaller clusters as well as on single, high-symmetric, larger ones. Alternatively, unbi- 
ased structure optimizations for larger clusters are all based on more approximate total- 
energy methods like the EAM or related methods , Si 9 ! 10 ! 11 ! 12 ! 13 ! 14 ? 15 ! 16 ! 17 whereby only for 
iV < 55 the structure has been optimized completely unbiased, or the very simple n-body 
Gupta , ^i 1 ?! 2 ?! 21 ! 22 Sutton-Chen,-^ 3 - Finnis- Sinclair,^ Murrell-Mottram,^ or Morse potential 
has been used, whereby unbiased structure optimizations up to N = 80 have been car- 
ried through.— Only for very simple and material-unspecific potentials, like the Lennard- 
Jones potential, largely unbiased structure optimizations up to N ~ 150 have been car- 
ried through. 24 Finally, Nijv clusters have also been studied with a tight-binding method,— 
whereby the structure only for clusters with N < 10 was optimized unbiased. 

Of the experimental studies on Ni^ clusters , 2 ?i 3 ?i 31 i? 2 i 33 i? 4 i 3 ^i?^i 37 the chemical-probe 
experiment o 29l3 ° l31132133134135136 have given very valuable information on the structure of the 
smaller clusters. On the other hand, results for medium and large Ni clusters,— obtained 
by performing near-threshold photoionization and time-of-flight mass spectroscopy can be 
used in identifying particularly stable clusters and, subsequently, in providing information 
on growth modes. 

/,From the earlier studies on nickel clusters it has been found that highly stable clusters 
occur for N = 13 and iV = 55 for which the structures are multilayer icosahedra, i.e., 
the first and second Mackay icosahedra. Moreover, in several theoretical studie d 17 ' 22 ' 23 a 
fee truncated octahedron was found for Ni 38 which is in accord with experimental results. 
Thus, it can be assumed that for these three cluster sizes a 'structural consensus' has been 
achieved, which in turn provides a useful check for any further study. On the other hand, 
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for N = 75 some theoretical studie o 22 ' 23 have predicted Marks decahedron to be the most 
stable structure, in disagreement with experiment. Moreover, the cluster growth pattern for 
intermediate sizes, i.e., for 13 < iV < 55 and for iV > 55, remains unclear and the obtained 
results seem to depend strongly on the used approach. 

In the present study the structure and the energetics of the four most stable isomers for 
small and medium-sized Ni^ clusters with N from 2 to 150 have been determined for each 
cluster size using a combination of the EAM (for calculation of the total energy for a given 
structure), the variable metric/quasi- Newton method (for the determination of local total- 
energy minima), and our own Aufbau/Abbau method (for the determination of the global 
total-energy minimum). Except for two cases (N = 75 and N = 104; here we included the 
decahedron as a starting configuration) our structure-determination is completely unbiased. 
In particular, as we shall see, it is important to go beyond N = 80 (the upper limit in the 
earlier unbiased structure-optimizations) in order to identify growth patterns. In addition, 
we shall present results for clusters that are obtained as relaxed spherical parts of the fee 
crystal with N up to some 1000s. 

The paper is organized as follows. In Sec. II we briefly outline the embedded-atom 
method, and in Sec. Ill we present our structural-determination methods. The main results 
are given in Sec. IV, and a brief summary is offered in Sec. V. Finally, for the sake of 
completeness we mention that some preliminary results on clusters with up to N = 100 and 
results for the energetically lowest isomer were presented earlie r 3Sl39 and in the discussion of 
the results we shall use eV and A as units of energy and length, respectively. 

II. EMBEDDED- ATOM METHOD 

We use the EAM of Daw, Baskes, and Foiles^* 4 ^ 4 ^ 3 " for the calculation of the total energy 
of a given cluster with a given structure. The main idea of Daw and Baske o 4 ^ 41 is to split the 
total energy of a (metallic) system of interest into two components. The embedding energy 
is obtained by considering each atom as an impurity embedded into a host provided by the 
rest of the atoms. The remaining part is written as a sum of pair potentials. Accordingly, 

£tot = £F<(rf) + i £ o,, (/«>,,). (l) 
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where p\ is the local electron density at site i, Fi is the embedding energy, i.e., the energy 
required to embed an atom into this density, and <pij is a short-range potential between 
atoms i and j separated by the distance R^, 

•■A 7,.(-R..\ 

(2) 



_ Zi(Rij) Zj(R i j / 



Here the effective charges Zi(Rij) and Zj(Rij) depend on R^. The local electron density at 
site i is assumed being a superposition of atomic electron densities 

pi = E pliRi^ (3) 

where p^(Rij) is the spherically averaged atomic electron density provided by atom j at the 
distance R^. 

In the presen t appr oach, the atomic densities have been taken from Hartree-Fock calcu- 
lations of Refs. 44M5] . Moreover, adjustable parameters that define Fi and 0^- have been 
obtained by fitting to known bulk properties such as sublimation energy, lattice constant, 
the heat of solution of binary alloys and, additionally, to the universal equation of Rose,— 
which describes the sublimation energy of the most metals as a function of lattice constant. 
The values for pf, Fi and Zi that we have used are available in numerical form for Ni, Cu, 
Ag, Au, Pd and Pt. 47 

Our main reason for choosing the EAM is that the EAM provides a computationally 
efficient, parameterized many-particle method that allows for many calculations of also larger 
systems. Thus, with this method we can perform unbiased structure optimizations for 
clusters with well above 100 atoms, which is not possible using more accurate, parameter- 
free methods. Moreover, by comparing (see Sec. IV A) with results of more accurate studies 
for the smallest clusters we can estimate the accuracy of the method. Since the method 
has been developed first of all for macroscopic metallic systems, in particular the smallest 
clusters may be those for which the largest inaccuracies show up. 



III. OPTIMIZATION OF THE STRUCTURE 

Using expression (JTJ) we can calculate the total energy of any cluster with any structure as 
a function of atomic coordinates {Ri}, E tot (Ri, R 2 , . . . , Rn)- In order to obtain the closest 
local total-energy minimum we use the variable metric/quasi- Newton method.— We found 
that this was significantly more efficient than the conjugate gradient method. 
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For searching the global minima we have developed our own Aufbau/Abbau method. It 
consists of the following steps: 

1) We consider two cluster sizes with N and N + K atoms with K ~ 5 — 10. For 
each of those we study a set of randomly generated structures, iV ran ~ 1000. Using the 
quasi-Newton method the iV ran relaxed structures are identified and the structures of the 
lowest total energy selected. Each of the iV ran starting structures for a cluster with M atoms 
was generated using a random-number generator for positions within a sphere or a cube of 
volume V c i = (p ■ b nn ) 3 M, where b nn = 2.49 A is the nearest-neighbor distance of bulk Ni 
and p = 0.8, 1.0, 1.2, i.e. we considered slightly compressed, normal, and slightly expanded 
structures. We included the constraints that the smallest allowed inter-atomic distance was 
0.56 nn and each atom has to interact with at least two others. 

2) One by one, each of the M atoms is displaced randomly, and the closest local minima 
is determined. If the new structure has a lower total energy than the original one, this new 
one is kept, and the old one discarded. This is repeated approximately 500 — 1000 times 
(depending on cluster size). 

3) This leaves us with two 'source' clusters, Ni^ and Nijv+K with their lowest total 
energies. One by one an atom is added at a random position to the structure with N atoms 
(many hundred times for each size), and the structures are relaxed. In parallel, one by one an 
atom is removed from the structure with N + K atoms — for each intermediate cluster with 
N' atoms we consider all N' + 1 possible configurations, that one can obtain by removing 
one atom from the NiAr/+i cluster. /,From the two series of structures for N < M < N + K 
those structures of the lowest energies are chosen and these are used as seeds for a new set of 
calculations. First, when no lower total energies are found, it is assumed that the structures 
of the global-total-energy minima have been identified, and we proceed to larger clusters. 

Moreover, by keeping information on not only the single energetically lowest isomer, but 
more low-lying ones, we have been able to identify the energetically four lowest-lying isomers. 

Our method combines randomness with regularity and is significantly more efficient com- 
pared to the case when only random starting geometries are used. In addition, it is com- 
pletely unbiased. However, as for any structure-determination method, there is no guarantee 
that we have identified the true global total-energy minima, although we believe that we 
in most cases have very good candidates. But, unfortunately, in two cases, N = 75 and 
N = 104 (both cases correspond to decahedral structures, which often have been difficult to 
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find in unbiased structure optimizations), our approach failed and we had to include that in 
our study explicitly, and in another, N = 38, we only got the previously obtained structure 
after considerably more attempts than otherwise used. We add that similar problems have 
been observed and analysed in studies on Lennard- Jones clusters.— 

IV. RESULTS AND DISCUSSIONS 

A. Small clusters — accessing the accuracy 

The smallest clusters, i.e., Ni^ for 2 < N < 13, provide an excellent test for our approach 
of two reasons: The clusters are so small that the number of (meta-)stable structures is 
small, which in turn means that the global total-minimum within a given approach can be 
identified. And since the EAM has its foundations in considerations for infinite, extended 
metallic systems with largely delocalized electrons, the small clusters should be those that 
with most difficulties can be treated by our approach. Therefore, we have in Table|l]tabulated 
our obtained structural parameters for 2 < N < 8 and iV = 13 in comparison with results 
of other theoretical or experimental studies. 

The table shows that the structure, here given through the average bond length (r), is 
most often given within a couple of % of the results of experiment or of parameter-free 
electronic-structure calculations. Moreover, none of the other parameterized methods, i.e., 
the tight-binding approach or model potentials of the Sutton-Chen, Gupta, or effective- 
medium type, gives results that are systematically more accurate than our approach does. 
This is important since for the larger clusters that we shall discuss here it is not possible 
to apply accurate parameter-free methods when simultaneously attempting to optimize the 
structure. 

The table also shows that the relative ordering of the different isomers may differ for 
different methods. Hence, for N = 7 most parameterized methods predict the isomer 
to have a lower total energy than the structure, whereas experiment and parameter-free 
theoretical studies predict a reversed ordering. Therefore, it can be important to calculate 
more different, energetically low- lying isomers when studying theoretically the properties of 
not very small clusters — identifying a different isomer as the one of the lowest total energy 
may suggest a significant disagreement between different studies. 
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In order to emphasize this point further we show in Fig. ^ the structures of the two 
energetically lowest isomers for 6 < iV < 13. It is very clearly seen that in some cases the 
two isomers are indeed very different and identifying the 'wrong' one as the most stable one 
may give unexpected results. 

It may also surprise that for the two isomers for N = 7 that look very different we find 
a total-energy difference of only 0.03 eV/atom. To our knowledge there exist two DFT 
calculations of Nij cluster.^ Both studies predict for the two energetically lowest isomers 
a capped octahedron and a pentagonal bipyramid, respectively. The corresponding average 
binding energy per atom are found to differ by only 0.05-0.07 eV/atom, i.e., very similar to 
our findings. 

Table and Fig. Q give one further result for N = 7: Ni 7 is the smallest cluster that 
incorporates the fivefold symmetry axis, an important element of icosahedral and decahedral 
symmetries. 

For N = 8 our two energetically lowest isomers differ in energy by 0.016 eV/atom, in 
good agreement with the value of 0.02 eV/atom found by Desmarais et al. in a DFT 
study. 6: There are some fundamental differences between different ab initio studies on Nig 
cluster . 3 i 4 i 6 i 7 According to one study 4 - a cube is the most stable structure for Nig, whereas 
in other studies^ a cube was only the 5th isomer. Further, in another work£ the second 
isomer was found to be a bicapped trigonal prism (C 2v ) and a capped pentagonal bipyramid 
was only isomer number four. Michelini et a/.— has found for the eight-atom nickel cluster 
only one isomer with the tetrahedral symmetry T^. This GGA result is inconsistent with 
all other studies. It seems that the LSDA calculations of Desmarais et al— have the best 
agreement with both experimental and semiempirical studies. 

To our best knowledge there are only DFT studies?' 4 on Nijv, 9 < iV < 13, for iV = 13. 
At iV = 13 the first Mackay icosahedron (I/J is constructed. The second isomer is 1.16 eV 
less stable. This is a rather large value and the icosahedron as the most stable structure 
for Nii3 was found almost in all theoretical studies (the only exception being an extended- 
Hiickel study 52 ). Once again we obtain good structural agreement with ab initio studies.— ^ 
The obtained bond lengths of 2.36 (center to vertex) and 2.48 A (vertex to vertex) agree 
well with those derived by Calleja et al.fi i-e., 2.41±0.03 and 2.53±0.03 A, and with those 
of Reuse and Khanna,— , i.e., 2.23 and 2.34 A, respectively. Furthermore, when comparing 
the differences in the binding energies per atom between the Ni 13 and Ni 2 clusters, our value 
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of 1.525 eV/atom agrees very well with the corresponding value of 1.558 eV/atom obtained 
using the self-consistent DFT program SIESTA.— 

Closing this subsection we emphasize that for our purpose, the calculation of energetic and 
structural properties of Nijy clusters, the present EAM approach appears to be sufficiently 
accurate. 

B. Energetic properties 

In Fig. |21 we show the binding energy per atom for the optimized clusters for the four 
energetically lowest isomers together with the binding energy per atom for clusters that 
were constructed as relaxed, spherical parts of the fee crystal structure (with an atom at the 
center of the sphere) for N up to almost 2600. With few expectations the four total-energy 
curves are so close that their difference can not be resolved in the figure. On the other hand, 
except for N = 79, the /cc-derived clusters lie all below the other structures, indicating that 
for iV up to around 150 the fee structure is not playing a significant role as a structural motif 
(an exception occurs for N = 38 for a fragment of the fee crystal with the center not at an 
atom — this case is not shown in the figure). Fig. Efb) shows that the apparent saturation 
is not real: first when considering considerably larger clusters, the convergence towards an 
average binding energy around 4.45 eV/atom can be recognized. 

In order to identify the particularly stable clusters we shall consider two criteria. First, 
a cluster may be considered very stable if its binding energy per atom is much larger than 
that of the two neighboring clusters. This can be quantified through the stability function, 
E tot (N + 1.1) + E tot (N - 1.1) - 2E tat (N.l), where E tot (N.k) is the total energy of the ener- 
getically fc-lowest isomer of the Nijy cluster. This function, that has maxima for particularly 
stable clusters, is shown in Fig- El Here we can identify a large number of particularly stable 
clusters, i.e., so-called magic clusters. These are found for N = 4, 6, 10, 13, 15, 19, 23, 26, 
28, 32, 36, 39, 43, 46, 49, 55, 58, 61, 64, 67, 71, 75, 77, 79, 83, 86, 89, 92, 95, 99, 101, 104, 
108, 110, 114, 116, 119, 122, 125, 127, 129, 131, 137, and 147. The most pronounced peaks 
occur at iV = 13, 19, 23, 46, 49, 55, 71, 75, 116, 131, 137, and 147 and most of them contain 
an odd number of atoms and, as seen in Table m these clusters possess high symmetry, 
including the highest icosahedral symmetry for Ni 13 , Ni 55 , and Ni^y. 

According to our other criterion for a particularly stable cluster, such a cluster occurs if 
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the energy difference between the two energetically lowest isomers E to t(N.2) — E tot (N.l) is 
large. This energy difference together with the energy differences to the energetically higher 
ones are shown in Fig. 0] It is striking that the energetically higher-lying isomers become 
energetically less and less separated and, as seen when comparing with Fig. |2l many of the 
clusters that are particularly stable according to the first criterion is it also according to the 
second one. 

C. Structural properties 

A central issue of our work is to identify general properties of the Ni^ clusters from 
the large amount of information that we have obtained from the calculations. Instead of 
discussing the individual clusters, in particular their structure (a brief summary of the 
structures is given in the Appendix), we shall introduce different quantities that are devised 
to reduce the available information to some few key numbers. 

First we shall consider the overall shape of the clusters. As we discussed in our earlier 
report, 39 it is useful to examine the 3x3 matrix containing the elements 

1 N 

ht = — ^{Rn,s — Ro,s){Rn,t ~ Ro,t) (4) 
U l n=l 

with ui = lA being our length unit, and s and t being x, y, and z, and with 

1 N - 

R0 = ^T, R n (5) 

JV n=l 

being the center of the cluster. The three eigenvalues of this matrix, I aa , can be used 
in separating the clusters into being overall spherical (all eigenvalues are identical), more 
cigar-like shaped (one eigenvalue is large, the other two are small), or more lens-shaped (two 
large and one small eigenvalue). Moreover, the average of the three eigenvalues, (I aa ), is a 
measure of the overall extension of the cluster. Finally, for a homogeneous sphere with N 
atoms, the eigenvalues scale like iV 5//3 . Therefore, we show in Fig. El quantities related to 
these eigenvalues but scaled by iV~ 5 / 3 . 

In the figure we see that only some few clusters have an overall spherical shape (these are 
found for the energetically lowest isomer for N = 4, 6, 13, 26, 28, 38, 55, and 147 and for 
the next one for N = 54, 79, and 92), which all correspond to high-symmetry isomers (cf. 
Table |HJ) and, for the lowest-energy isomer, also to the class of magic clusters. Moreover, it 
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is interesting to notice that the average value follows more or less the same curve for all four 
isomers, i.e., the more compact ones (with relative low average value) are found for N just 
above 50 and above 140). Also the largest differences follow more or less the same curve, 
except for some few cases mainly for N below 40. Therefore, except when the eigenvalues 
all are very similar (i.e., the largest difference is very small, which occurs for iV around 50, 
100, and 140), the overall shape (i.e., lens- or cigar-like) is the same for all four isomers. 

In our earlier work 39 we showed that the structural development of the isomer of the lowest 
total energy with advantage could be monitored through so-called similarity functions. Thus, 
for a given cluster Nijv we define a radial distance for each of the atoms 

r n = \R n - R \ (6) 

and sort these in increasing order. 

Simultaneously we consider a large fragment of a fee crystal or, alternatively, a large 
cluster of icosahedral symmetry. Also for these we define a radial distance for each atom, 
r' n , that also are sorted. In order to compare the given cluster with any of those systems we 
calculate subsequently 
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1 N -,1/2 
iV n=l 



(7) 



giving the similarity function 



S = 1 1 / ( 8 ) 



(ui = lA), which approaches 1 if the Ni^r cluster is very similar to the reference system. In 
Fig. we show the resulting functions in four cases, i.e., when comparing with the relaxed 
Ni 30 g cluster of icosahedral symmetry, and when comparing with three fragments of the 
fee crystal differing in the position of the center (the position of an atom, the middle of a 
nearest-neighbor bond, and the center of the cube, respectively). 

We see that both clusters that clearly resemble fee fragments and such that resemble 
icosahedral clusters can be identified. Of obvious reasons, the latter is the case for N around 
13, 55, and 147, but it is interesting to notice that these clusters are not singular: clusters 
with a large similarity to icosahedral clusters are found over a larger range of clusters (~10) 
around the truly icosahedral cluster. Similarly, clusters that are close to being fragments of 
the fee crystal can be identified around N = 13, 28, 38, and 79. Moreover, it is surprising 
to observe that all four isomers for a given size show very similar similarity functions. 
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Also the radial distances themselves can give some interesting information. These are 
shown in Fig. [7| for the four different isomers as functions of cluster size. Up to iV around 
40, only few trends can be identified (with the exception of iV = 13), but for N just above 
50 a clear tendency towards shell-construction can be seen. Whereas this transition towards 
shell-construction is fairly smooth as a function of cluster size, independent of whether we 
start from the larger or the smaller clusters, in the range 75 < N < 80 abrupt changes 
occur also towards some shell-construction. That this size range is special shall be discussed 
in more detail below. Also for N close to 115 this tendency is seen, but most dramatic 
are the results for clusters around 140 and above where, for all four isomers, a very clear 
shell-construction is identified. Comparing with Fig. |U] we see that this is correlated to the 
construction of icosahedral-like structures. Thus, since the infinite crystal possesses the fee 
structure, the clear shell-construction in this size range will have to disappear again for 
larger clusters. 

In total, it can not surprise that when analysing the structures of the individual clusters 
in detail we can identify many structural motifs that are characteristic for clusters of par- 
ticularly high stability also for clusters with somewhat different sizes. Moreover, such an 
analysis^ shows also that many of our structures for the smaller N (these are the only ones 
that have been the subject of earlier unbiased structure-optimizations) are close to those 
that have been found in earlier studies. Here, however, we shall largely omit the discussion 
of the individual structures but instead explore the more general trends. 

In Fig. El we show some few quantities related to the structure of the clusters, i.e., coor- 
dination numbers and average bond lengths. Here, we define two atoms as being bonded if 
their interatomic distance is less than 3.00 A, which is the average of the nearest-neighbor 
distance (2.49 A) and the next-nearest-neighbor distance (3.52 A) in bulk Ni. Moreover, 
we distinguish between inner atoms with a coordination number of 12 or larger and surface 
atoms with a coordination number less than 12. 

Fig. lEfa) shows the average coordination number as function of N, where an obvious 
increase as function of size is observed, but with additional oscillations in particular for the 
smallest clusters. This function shows clearly a saturation towards the bulk limit of 12, 
although one has to remember that even for the largest cluster of our study 94 out of 150 
atoms are characterized as surface atoms. 

Fig. lEfb) shows the minimum atomic coordination for each cluster size from the range 
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N = 2 — 150, which is useful information for chemical-probe experiments due to the enhanced 
chemical activity of clusters containing low-coordination atoms. In addition, the existence 
of low-coordination atoms (i.e., coordinations of 3 or 4) could point to the occurrence of 
a cluster growth, where extra atoms are added to the surface of the cluster, whereas the 
higher minimum- atomic-coordination numbers (e.g. 5 or 6) could indicate a growth where 
atoms are inserted inside the cluster, or, alternatively, upon a strong rearrangement of 
the surface atoms. These considerations follow the so-called "maximizing the minimum 
coordination" . 16 i 17 

Fig. |HJ(c) shows the average bond length as function of cluster size. The dashed line 
corresponds to the bulk value of 2.49 A. The average bond length for all the structures with 
the exception of Ni 31 , Ni 35 , Ni 71 , and Ni 86 is smaller than the bulk value. Such a reduction 
is in accordance with the results of most other theoretical studies, but it is interesting to 
observe that this property approaches the bulk value faster than the other quantities of the 
figure. 

Important experimental work in the determination of the structure of individual Nijv 
clusters has been presented by Riley et a / . 30 i 31 i 32 i 33 i 34 i 35 using chemical-probe experiments. 
Thereby, the number of N 2 molecules that is bonded to a given cluster is determined, and in 
order to relate this number to the structure of the cluster, the following empirical rules have 
been formulated : 30 ! 31 ' 32 (i) nickel atoms with a coordination of four and less will bind two 
nitrogen molecules at saturation, (ii) nickel atoms with coordination of five to eight will bind 
only one nitrogen molecule at saturation, (Hi) nickel atoms with a coordination of nine will 
bind one nitrogen molecule weakly or not at all (this kind of binding has not been observed 3 ^ 
for N > 49), and (iv) nickel atoms with a coordination more than nine do not bind any 
nitrogen molecules. These rules allow us in turn to estimate the number of N 2 molecules 
that will bond to a given cluster. We stress that the chemical-probe experiments provide 
information on the complete system, i.e., cluster plus adsorbants, whereas we consider the 
naked clusters. Although this difference may result in some deviations between theory and 
experiment, a comparison should give useful information on the clusters. This is supported 
by the studies of Parks et a/ . 30 ' 31 that show that in many cases the structures of naked and 
N 2 -covered clusters are very similar. 

Doing so, we obtain the results of Fig. El where we also compare with the results of other 
theoretical works. The experimental work covers the range 3 < N < 72, whereas the former 
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theoretical works cover 4 < iV < 80 (Ref. H), 3 < N < 23 (Ref. 2^), and 24 < iV < 55 
(Ref. [l7[), obtained using the Sutton-Chen potential,— the Finnis-Sinclair potential and 
the effective- medium theory.— Fig. Efa), where we show the results for the whole range of 
our study, i.e., 2 < N < 150, shows an interesting aspect that obviously is beyond most 
of the previous experimental and theoretical studies, i.e., that up to around iV ~ 55 the 
number of binding sites increases more rapidly with TV than for N above ~55, i.e., the curve 
is roughly two piecewise linear curves. The experimental studies [cf. Figs. Efd) and (e)] do 
in fact show a small tendency in this direction, which, on the other hand, is not found in the 
earlier theoretical work [cf. Fig. EUc)]. Figs. Efb)-(e) show a comparison of the experimental 
results with those of the four different theoretical studies to this issue. It is clear that our 
approach is at least as accurate as any of the other theoretical ones. Moreover, the Sutton- 
Chen potential seems to predict a larger scatter of the number of binding sites as a function 
of size compared to any other approach. 

Whereas icosahedral motifs are found to play an important role for very many values 
of N, we found that for iV around 75 also other structural motifs are important. The 
chemical-probe results for N2 uptake, both the experimental ones (notice, however, that the 
experimental studies are limited to iV < 72) and the simulated, theoretical ones, suggest 
also that in this range of N values more structural motifs compete: for these values of N the 
results of Fig. El are particularly scattered. Moreover, as we shall see in the next subsection, 
also growth processes suggest that these N values are special. Therefore, we shall study the 
structural properties of those clusters in some more detail. 

For cluster sizes ranging between 70 and 80 atoms two more closed-shell structures other 
than icosahedral ones can be identified, i.e., the (2,2,2) Marks decahedron (giving a 75-atom 
cluster) and a 79-atom fee cluster, which is the truncated centered octahedron. Due to their 
quasi-spherical shape they have a particularly low surface energy. Moreover, the cluster 
sizes between 70 and 80 are very different from the values of closed-shell icosahedral clusters 
that have 55 or 147 atoms. Accordingly, one may expect a strong competition between 
these three main morphologies in this size range. Our calculations support this suggestion. 
Fig. El shows the lowest total energies for the decahedral and fee clusters relative to the 
icosahedral ones for 72 < N < 82. The Marks decahedron and its derivatives dominate over 
the icosahedral structures for 73 < N < 79. At N = 79 both decahedral and fee motifs 
are more stable than the icosahedral one. The decahedron at N = 79 is only 0.04 eV more 
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stable than the truncated octahedron. Finally, three of four lowest-energy isomers at iV = 75 
are decahedra and for N = 76 — 78 there are only decahedral clusters among the four first 
isomers. 

This result suggests that a fundamental structural change takes place for iV around 
75. This, however, is only partly the case: for larger N values the icosahedral packing is 
again stabilized (an exception occurs for N around 104, where decahedral structures are 
found). In order to verify this we computed the total energies for several optimized closed 
fee and decahedral structures, i.e., truncated noncentered octahedra for N = 116 and 140, a 
noncentered octahedron and a (3,2,2) Marks decahedron for N = 146, and a cuboctahedron 
and an Ino decahedron for N = 147. In all cases the icosahedral structures were found to 
have lower total energies. The fee clusters for N = 116, 140 and 146 and m-Dh 146 were 0.85, 
2.01, 5.64 and 2.29 eV, respectively, less stable than the icosahedral-like structures. The 
three-shell Mackay icosahedron was 5.69 and 4.71 eV more stable than the cuboctahedron 
and the Ino decahedron, accordingly. Thus, our calculations point strongly to icosahedral 
growth up to cluster sizes of at least 150. This statement is supported by mass-spectroscopy 
studies of nickel clusters^ which suggest icosahedral growth over a wide cluster-size region 
from 50 to 800. 

D. Growth patterns 

One of the central issues is to identify how clusters grow, i.e., whether some fundamental 
principles determine how (and whether) the cluster with N atoms is derived from the one 
with N — 1 atoms. We shall here study this issue from a completely static point of view, i.e., 
dynamic and kinetic effects, that most likely are important in experiment, will be neglected 
here. Basic information about this issue is given by the structures of the magic clusters. The 
majority of the structures of these are icosahedral. Thus our study supports predominantly 
icosahedral growth of nickel clusters. There are two different ways of icosahedral growth. 
Thus, the first Mackay icosahedron, Ih 13 , contains 20 triangular faces joined by 30 edges 
and 12 vertices. In the first case, MIC/Mackay or fee covering , 11 i 55 i 56 i 57 i 58 additional surface 
atoms are placed on top of the edges and the vertices. This kind of growth results ultimately 
in the second Mackay icosahedron, II155. Another possibility, the TlC/polyicosahedra l 111 ^ 
or hep or face-capping growth (see Ref. and references therein), corresponds to adding 
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atoms on top of the atoms at the center of each face (T sites) which leads, after complete 
covering, to the 45-atom rhombic tricontahedron. The TIC structures are the structures 
with the shorter average bond length, but at the same time they possess the larger strain 
energy. They are favored at the beginning of a covering until some critical size (depending 
on details of the atom-atom interactio n 57 ! 58 ) after which the MIC structures become lower 
in energy. 

To explain the icosahedral shell filling in Mg and Ca clusters Martin et a/.— proposed 
the so-called umbrella model. According to this model enhanced stability is achieved each 
time a vertex and its five surrounding faces are covered (capping umbrellas). The capping 
umbrellas for MIC and TIC structures lead to different magic numbers. By covering of Ih 13 
the MIC umbrellas are completed at iV = 19, 24, 28, 32, 36, 39, 43, 46, 49 and 55, whereas 
the TIC umbrellas give magic numbers for 19, 23, 26, 29, 32 and 34. Comparing with our 
set of magic numbers (see above), we find that growth from Ni 13 to Ni 2 7 is icosahedral and it 
can be described according to the umbrella model. For N > 27 our magic numbers coincide 
with the MIC magic figures. This confirms seemingly the suggestion^ that at iV = 27 there 
is a transition from TIC to MIC growth. However, a detailed analysis of our structures 
shows that our structures for Ni 281 and Ni 32 .i are not icosahedral. Instead, the most stable 
icosahedral structure appears first at N = 35 whereas the structures according to the MIC 
growth or to fee packing is found up to N = 56 with a single exception for N = 38. Thus the 
transition from anti-Mackay to Mackay covering occurs in the region of cluster sizes from 
28 to 34. Cluster growth in this range is complicated and there is a competition between 
icosahedral, polyicosahedral and tetrahedral structure. Probably this non-trivial growth was 
responsible for the difficulties in interpreting results of nitrogen-uptake experiment o 3 ° l31132 
for 29 < N < 48. 

At N = 55 we have the next complete icosahedron. According to the TIC growth model 
we should then have magic number at iV =58, 61, 64, 67, and 71, 11 which is in excellent 
agreement with our results. At N = 71 the TIC umbrella growth is completed. The MIC 
umbrella growth is completed at iV = 71, 83, 92, 101, 110, 116, 125, 131 and 137, again in 
agreement with our calculated magic numbers. 

Our study allows a more detailed analysis of possible growth processes that ultimately 
will confirm the complexity of the growth. Thus, in order to quantify whether the cluster of 
iV atoms is related to that of N — 1 atoms we consider first the structure of the lowest total 
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energy for the (N — l)-cluster. For this we calculate and sort all interatomic distances, di, 
i = 1, 2, ■ • • , N ( N 2 ~ 1 > . Subsequently we consider the iV fragments of the A^-cluster that can 
be obtained by removing one of the atoms and keeping the rest at their positions. For each 
of those we also calculate and sort all interatomic distances d[, and subsequently calculate 



q 



2 N(N- l)/2 i 2 



(9) 



Of the A" different values of q we choose the smallest one, q min and calculate the similarity 
function 

S = - 1 —— (10) 

1 + Qmin /Ul 

(ui = lA) which approaches 1 if the Nijv cluster is very similar to the Nijv_i cluster plus an 
extra atom. 

This function is shown in Fig. HIT a). We see indeed that for iV up to around 40, S is 
significantly different from 1, giving support for the consensus that in this range the growth 
is complicated. For larger A" we see again that around iV ~ 75 some larger structural changes 
take place. Fig. |3] shows, however, that the energetically lowest isomers differ only little in 
energy. Therefore, one may ask whether isomers with a slightly higher energy play a role 
in the growth process. Accordingly, we examined whether the structure of the energetically 
lowest isomer of Ni^ would resemble any of the four energetically lowest isomers for Nijv-i in 
the same sense as above, i.e., we use the quantity q in quantifying the structural difference. 
Out of the four different values for q (for the four different isomers) we choose the smallest 
one and construct the similarity function from this. The resulting function is shown in 
Fig. fTTT b). and compared with Fig. HTT a) it is clearly seen that in particularly in the region 
N > 50 the similarity function now approaches 1, except for the clusters around A" = 70 — 75, 
i.e., the region where we already in analysing other properties have found atypical behaviors. 
In Fig. lTTT c) we show which of the four energetically lowest isomers of Ni^r-i leads to largest 
similarity with the Nijv cluster, and from the fact that this only for the smaller A" most often 
is the first isomer, we learn that the growth process is a complicated process where more 
different isomers are important, i.e., not only the energetically most stable one, although 
these may resemble each other (cf. Fig.|5J). This means that it is difficult to imagine that the 
growth occurs as a one-by-one addition of atoms to the energetically most stable isomers. 
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V. CONCLUSIONS 



In summary, we have determined the four energetically lowest isomers of nickel clusters 
Nijv in the range 2 < N < 150, using a combination of the embedded-atom method in 
the version of Daw, Baskes and Foiles (for the calculation of the total energy for a given 
structure), the variable metric /quasi- Newton method (for the determination of the closest 
total-energy minimum), and our own Aufbau/Abbau method (for the determination of the 
global total-energy minimum). Although the calculations provide structural information 
for each individual cluster, separately-^ that may be useful information for experimental 
studies, we have refrained from discussing these separately (some information is, however, 
given in the Appendix), but instead focused on identifying general trends for Nijv clusters, 
and reducing the information on the individual clusters to some few key quantities like 
total energy (per atom), overall symmetry and shape, average bond length and coordination 
number, and similarity with other structures (fee fragments, icosahedral clusters, or smaller 
clusters) . 

Since our total-energy approach (the embedded-atom method) has its roots in arguments 
for extended solids with delocalized atoms, we first verified that it could be applied also 
for the fairly open structures of our interest. We found that for the absolutely most 'crit- 
ical' cases, i.e., the smallest clusters, our results are in excellent agreement with available 
theoretical and experimental information making us believe that our approach is accurate. 
Furthermore, we have also tried to use the parameterization of Voter and Che n 61 ' 62 of the 
EAM, which has been produced in order to give accurate descriptions of also finite and more 
open systems, for some of the smallest clusters. To our satisfaction, no structural changes 
were found and, therefore, we did not pursue this approach. 

Our study predicts a number of particularly stable clusters, i.e., magic numbers, which 
in many cases are in agreement with the prediction of other studies and in some cases with 
'common expectations' although we stress that in our study they were found in a completely 
unbiased approach. These magic numbers were clearly visible both in the 'stability function' 
and in the total-energy difference between the energetically lowest and higher-lying isomers, 
whereas for the energetically higher-lying isomers, the differences in the total energy become 
small. 

We found also that even for large values of N (above 2500) the total energy per atom 
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has still not converged to the bulk limit. Similarly, the average coordination number is for 
clusters with N around 150 still far from the bulk limit, whereas the average bond length 
has come close to the bulk limit for this size. 

In order to study the structural properties we analysed the eigenvalues of the matrix 
containing the moments of inertia and introduced so-called similarity functions. The latter 
indicated that all the four isomers we have studied here are very similar in most cases, and 
that roughly spherical clusters were found mainly for the energetically lowest isomer but in 
some cases also for the second-lowest one, and that these often correspond to particularly 
stable structures. With the similarity functions we could identify clusters with a predom- 
inantly icosahedral structure and such with an essentially /cc-like structure. Whereas the 
former occurred for iV around 13, 55, and 147, the latter was found, e.g., for iV around 79. 
It was interesting to see that these structural motifs were built up over a larger range of N 
values and were, accordingly, not limited to those singular values of N. 

By analysing the distribution of radial distances as a function of cluster size we could 
identify certain size regions, e.g., N around 50, 75, 115, and above 140, where shell structures 
were clearly recognizable. This property, as well as several others, showed, furthermore, that 
all four isomers for a given N have very similar structures. 

Chemical-probe experiments have been used in identifying the structure of the individual 
clusters, and for this we found an excellent agreement. Our study shows that for clusters 
with N around 50-60, a change occurs, so that above this size the number of N2 binding 
sites depends much weaker on N than below this size. This change is just on the boarder 
of the cluster sizes that have been studied experimentally and theoretically so far. 

We used the magic numbers in analysing different models for cluster-growth processes. 
Generally, the cluster growth is according to multilayered-icosahedral or layer-by-layer 
growth. We found that the one-, two-, and three-shell Mackay icosahedra are the most 
stable structures for 13-, 55-, and 147-atom nickel clusters (the latter for the first time 
in unbiased studies). This result agrees well with the suggestion in Ref. that Mackay 
icosahedra are the most stable structures up to iV = 2300. But the cluster-growth pattern 
between two closed-shell icosahedra is not simple. It is predominantly icosahedral with is- 
lands of fee growth for Ni 6 , Ni 38 , tetrahedral for Ni 2 s and decahedral for Ni 73 -N 7 9 and Ni 10 4. 
For 14 < iV < 27 the growth is polyicosahedral. Moreover, we found a growth according to 
/jcp/anti-Mackay for 57 < iV < 72 and for 148 < iV < 150. From iV = 35 to iV = 56 (with 
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the exception of N = 38) and from N = SO to N = 147 (except for N = 104) the MIC 
cluster growth was found. Thus, the change-over from the hep-type of cluster growth to the 
fee one always happens within a transition region of cluster sizes and not immediately. 

Moreover, whereas the growth process only with difficulties could be identified for the 
smaller N (due to the lack of similarity between the clusters), for larger N we found that not 
only the energetically lowest isomers could be important in describing the growth. Thus, 
for all ranges of N, the growth is non-trivial. 

A special case is that for N around 75. More different properties suggested that in this 
range special things occur. Analysing the total energy for different structures in this size 
range we found indeed that more structures were energetically close and, accordingly, would 
play a role for those clusters. Therefore, for these values of N there is a rapid and dramatic 
change in structure as a function of N. 
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APPENDIX A: STRUCTURE OF THE CLUSTERS 

In this appendix we shall present the structure of some of the clusters without entering a 
detailed discussion of each clusters individually. We shall use the notation N.k for the Nijv 
cluster of the fc-lowest total energy. We add that the nuclear coordinates and total energies 
for the optimized clusters can be obtained from the authors upon request. 

3.1: equilateral triangle. 4.1: tetrahedron. 5.1 trigonal bipyramid. 6.1 octahedron. 
6.2: pentagonal bipyramid without a single pentagonal vertex. 7.1: pentagonal bipyramid. 
7.2: capped octahedron. 8.1: bisdisphenoid. 8.2: capped pentagonal bipyramid. 9.1: 
bicapped pentagonal bipyramid. 9.2: tricapped trigonal prism. 10.1: tricapped pentagonal 
bipyramid. 10.2: trapezoidal antiprism. 11.1: quadricapped pentagonal bipyramid. 11.2: 
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two pentagonal bipyramids with a common face. 12.1: icosahedron without one atom at 
one vertex. 12.2: low symmetry. 12.3: trapezoidal bipyramid with capped trigonal faces. 
12.4: truncated trigonal bipyramid. 13.1: Mackay icosahedron. 13.2, 13.3, 13.4: similar 
structures as 13.1, but with the move of a vertex atom to a face. 14.1: the structure of 13.1 
with the capping of one of the trigonal faces. 14.2: similar, but with the 14th atom on the 
top of a pentagonal edge. 14.3: bicapped hexagonal antiprism. 14.4: tricapped pentagonal 
prism. 15.1: bicapped icosahedron. 15.2: centered bicapped hexagonal antiprism. 15.3, 
15.4: capped pentagonal prism. 

16.1: tricapped icosahedron. 18.1: double icosahedra without one atom in one pentagonal 
ring. 18.2: double icosahedra without one atom in a capping vertex. 19.1: double icosahe- 
dron. 19.2, 19.3, 19.4: derivatives of the double icosahedron by moving of an atom from an 
apex or from a pentagonal ring to a face-capped or a bridging position. 20.1: double icosa- 
hedron with an atom added to the a pentagon in a bridging position. 21.1, 21.2, 21.3, 21.4: 
derivatives of the double icosahedron with two additional atoms at bridging or capping po- 
sitions. 22.1, 22.2: triple icosahedron without one atom either from a pentagonal ring (22.1) 
or from a capping position (22.2). 23.1: triple icosahedron. 24.1, 24.2: triple icosahedron 
with an additional atom. 25.1: three sequentially interpenetrating icosahedra whose vertices 
are located on one line. 26.1: six interpenetrating double icosahedra (a polyicosahedron). 
27.1: the polyicosahedron plus an extra atom at a pentagonal cap. 28.1: tetrahedron. 

In general, the structures for 29 < iV < 34 do not display any dominant morphology - 
there is a competition between tetrahedral, polyicosahedral, and icosahedral packings. One 
can consider this size range as a transition range from hep to fee growth. 

29.1: distorted tetrahedron. 30.1: a polyicosahedral morphology. 32.1: derived from the 
28.1 tetrahedron by adding four surface atoms. 35.1: MIC/Mackay structure. 36.1, 37.1: 
contain icosahedral motifs. 38.1: truncated octahedron. 38.2, 38.3: icosahedral fragments. 
39.1: the second Mackay icosahedron (55.1) without a 16-atom cap. 39.4: a fee cluster. 

49.1: the second Mackay icosahedron (55.1) without one face. 51.1, 52.1, 53.1, 54.1: 
obtained from 49.1 by adding atom by atom until 55.1: second Mackay icosahedron. 54.2: 
as 55.1 but without the central atom. 56.1: as 55.1 but with an extra atom at the center of 
a triangular face. N — 56 — 71: TIC/anti-Mackay growth in accordance with the umbrella 
model. 

75.1: (2,2,2) Marks decahedron. 79.2: a fee cluster. 104.1: related (2,3,2) Marks decahe- 
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dron. 147.1: three-layer icosahedron. 
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TABLE I: The symmetry (Sym.) and averaged bond length ((r)) for smaller Nijv clusters from dif- 
ferent experimental and theoretical studies. EAM, TBMD, FS, SC, AI, EMT, nG, and exp denotes 
Embedded-Atom calculations, tight-binding molecular-dynamics calculations, Finnis-Sinclair po- 
tential, Sutton-Chen potential, ab initio calculations, effective-medium-theory calculations, many- 
body Gupta potential, and experiment, respectively. For some values of N, more different meta- 
stable isomers have been found, and they have then been characterized by an additional number 



(e.g., 7.1, 7.2, 7.3, and 7.4, with 7.1 being the stabler one, 7.2 the second-most stable one, a.s.o.) 
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TABLE II: Point groups of the optimized clusters. 
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FIG. 2: Binding energy per atom as a function of size for the four energetically lowest isomers for 
N up to 150 (continuous curves) together with that of relaxed, spherical parts of the fee crystal 
structure with an atom at the center (closed circles). The horizontal dashed lines mark the value 
for bulk Ni as obtained with the same method. 
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FIG. 3: The stability function as a function of cluster size. 
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FIG. 4: The total-energy difference between the two energetically neighboring isomers as a function 
of cluster size, (a) shows the difference between the second and the first, (b) between the third 
and the second and (c) between fourth and the third isomer. 
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FIG. 5: Different properties related to the eigenvalues I aa . In the upper panel we show the average 
value together with points indicating whether clusters with overall spherical shape (lowest set of 
rows), overall cigar shape (middle set of rows), or overall lens shape (upper set of rows) are found 
for a certain size. Moreover, in each set of rows, the lowest row corresponds to the energetically- 
lowest isomer, the second one to the energetically second-lowest isomer, etc. In the lower panel we 
show the maximum difference of the eigenvalues for the four different isomers. 
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FIG. 6: Each panel shows the similarity function for all four isomers when comparing with (a) 
an icosahedral cluster, and (b-d) a spherical fragment of the fee crystal when the center of the 
fragment is placed at (b) the position of an atom, (c) the middle of a nearest-neighbor bond, and 
(d) the center of the cube, respectively. 
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FIG. 7: Each panel shows the radial distances (in A) as a function of cluster size, i.e., each small 
line represents (at least) one atom with that radial distance. The energy of the isomers increases 
from the top to the bottom. 
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FIG. 8: (a) the average coordination number, (b) the minimum coordination number, and (c) 
the average bond distances as functions of cluster size. The dashed lines in (a) and (c) show the 
corresponding bulk values for nickel. 
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FIG. 9: The (open circles) experimentally determined number of binding sites for Ni/v clusters 

nnn 

from Refs. 33.34.35] as function of N in comparison with theoretical calculated numbers (closed 
circles) from (a,b) the present work, (c) from 23], (d) from 25j], and (e) from jl 71 ] . 
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FIG. 10: The relative stability of the icosahedral, decahedral and fee structures as function of 
cluster size for 72 < N < 82. 
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FIG. 11: (a) and (b) show the similarity functions that describe whether the cluster with N atoms 
is similar to that of N — 1 atoms plus an extra atom, when (a) only considering the lowest-energy 
isomer for the (N — l)-atom cluster and (b) considering all four isomers for that cluster, (c) shows 
which isomer in the latter case is most similar to the one of TV atoms. 
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